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Break-up of excited ^^C in three equal-energy a particles 
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The fragmentation of quasi-projectiles from the nuclear reaction "'"Ca+^^C at 25 MeV/nucleon 
was used to produce excited states candidates to a-particle condensation. The methodology relies 
on high granularity A-k detection coupled to correlation function techniques. Under the assumption 
that the equality among the kinetic energies of the emitted a-particles and the emission simultaneity 
constitutes a reliable fingerprint of a condensation, we identify several tens of events corresponding 
to the deexcitation of the Hoyle state of ^^C which fulfill the condition. 

PACS numbers: 25.70.-z heavy ion reactions at low and intermediate energies, 
25.70.Pq fragmentation in nuclear reactions, 
67.85.Hj Bose-Einstein condensates, 
21.60.Gx Cluster model, nuclear structure. 



Bose-Einstcin condensation is known to occur in 
weakly and strongly interacting systems such as dilute 
gases and liquid ^He [ij. Since more than ten years it is 
also theoretically shown that for symmetric nuclear mat- 
ter, below a critical density, a-particle condensation is 
favored [il-Q- It is expected to occur at densities smaller 
than a fifth of the nuclear saturation density. At higher 
densities the 2-nucleon deuteron condensation prevails 
over the 4-nucleon condensation. This new possible phase 
of nuclear matter may have its counterpart in low-density 
states of self conjugate lighter nuclei, in the same way as 
superfluid nuclei are the finite-size counterpart of super- 
fiuid nuclear and neutron matter. This means that under 
some circumstances, the alpha condensation, i.e. bosonic 
properties, might dominate over the nucleon properties 
even in finite nuclei. Thus, by showing that the Hoyle 
state (i.e. the first excited state 0+ at 7.654 MeV of 
i^G) and the sixth 0+ state at 15.097 MeV of ^^O are 
described by a-particle condensate type functions, Refs. 
[1, 0] advance the idea that these states are candidates 
to observe a condensation. A common feature of these 
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states is their diluteness. For instance, calculations Q 
and recent experimental data 0) show that the rms ra- 
dius of the Hoyle state exceeds by 45% the radius of ^^C 
in its ground-state. On the other hand the Hoyle state is 
known to p lay a decisive role in stellar nucleosynthesis of 
^^G |8l4lC| with the assumption that the only contribu- 
tion to the alpha-particle width comes from the ^Be^.s. + 
a reaction |ll| . At present an experimental upper limit of 
4% exists for the contribution from the direct 3a channel 
to the alpha width [l^]. A phase space calculation sug- 
gests that the probability ratio of the 3a to the a-|-^Be 
decay process is 5 xlO~* [l3|. However such a factor 
may be enhanced in relation with the peculiar properties 
of the Hoyle state previously discussed. 

According to the present understanding, the compati- 
bility of a nuclear state with a-particle condensation may 
be judged upon its excitation energy close to the Na 
threshold, the emission simultaneity and both the low 
kinetic energy and kinetic energy dispersion. Then, it be- 
comes clear that probably the most appropriate method- 
ology should involve high velocity reaction products in 
the laboratory detected by a high granularity-high solid 
angle particle array, as needed in multifragmentation re- 
action studies [12] • 

The Letter reports on results obtained for the nuclear 
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reaction °Ca+^^C at 25 MeV per nucleon incident en- 
ergy performed at INFN, Laboratori Nazionali del Sud 
in Catania, Italy. The beam impinging on a thin carbon 
target (320/ig/cm^) was delivered by the Superconduct- 
ing Cyclotron and the charged reaction products were 
detected by the CHIMERA in multidetector [H]. The 
beam intensity was kept around 10'' ions/s to avoid pile- 
up events. CHIMERA consists of 1192 silicon-CsI(Tl) 
telescopes mounted on 35 rings covering 94% of the solid 
angle, with polar angle ranging from 1° to 176°. Among 
its most interesting characteristics, we mention the low 
detection and identification thresholds for light charged 
particles (LCP) and the very high granularity at forward 
angles. The mass and charge of the detected nuclei were 
determined by standard time of flight (TOE) for LCP 
stopped in silicon detectors and AE — E {Z > 5) and 
shape identification {Z < 5) techniques for charged prod- 
ucts stopped in CsI(Tl). In particular ^Be nuclei (two 
equal-ener gy a s hitting the same crystal) were identified 
in CsI(Tl) [1^. The energy of detected nuclei was mea- 
sured by the Si detectors calibrated using proton, carbon 
and oxygen beams at various energies ranging from 10 
to 100 McV. For Z ~ 2, dedicated energy calibrations 
of the fast component of CsI(Tl) light was realized us- 
ing the time of flight (TOE). Though high-quality TOE 
charts allowing for a direct calibration of the Csl light ex- 
ist only for 60% of the total number of modules, we have 
finally managed to calibrate more than 95% of modules 
from 1° to 62°. The modules for which TOE information 
was poor or missing were calibrated by comparing the 
fast component distribution with the benchmark distri- 
bution of the corresponding ring built out the telescopes 
with excellent TOE. 

Invariant Vpar — Vper plots show the binary character of 
the collisions with the formation of two emitting sources, 
a quasi-target (QT) and a quasi-projectile (QP). These 
latter may be easily separated according to the smaller 
or larger velocity of reaction products with respect to 
Vproj l'^- After a first event filtering according to the total 
a multiplicity {ma > 3) and the total detected charge 
{Ztot £ Zp + Zt = 26), we have focused exclusively on 
QP decay products with > 3. 

The power of correlation functions was extensively 
discussed in the literature [iB-El- They can furnish 
space-time information taking advantage of proximity ef- 
fects induced by Coulomb repulsion, reveal excited states 
(whose decay produces the detected correlated parti- 
cles) and enlighten any production of events or subevents 
with specific partitions. In this Letter various correla- 
tion functions have been used: two-particle correlation 
functions to first judge the quality of energy calibrations 
and multi-particle correlation functions to identify pos- 
sible a-particle condensation states and define their de- 
excitation characteristics. 

Two- or multiple-particle correlation functions, with 
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FIG. 1: (Color online) Correlation functions in total kinetic 
energy for a — a (top panel) and in excitation energy for Q-d 
(bottom panel) emission. Full and open symbols correspond 
to different ways of calculating the angle under which the 
detected particles were emitted (see text). The arrows in the 
insets point the corresponding centroids of the Breit-Wigner 
distributions (solid lines) fitted on the peaks. 

only one variable E are defined as: 

^ uncorr } 

where the role of the generic variable E is equivalently 
played by the total kinetic energy of the particles of inter- 
est in their center-of-mass frame Etot or by the excitation 
energy of their emitting source/state, E^x = Etot — Q- 
In Eq. (dl Ycorr, the correlated yield spectrum, is con- 
structed with the considered particles in the same event 
and Yuncorr stands for the uncorrelated yield spectrum 
constructed by taking particles in different events. In 
addition to these, by comparing the width of the peaks 
of the correlation function with the natural width of the 
excited states, one can estimate the distortions caused 
by the finite granularity. Fig. [1] presents the two-particle 
correlation functions for a — a (a) and a-d (b) in terms of 
total kinetic energy and, respectively, excitation energy. 
The error bars have been calculated taking into account 
the statistical errors on both the correlated and uncor- 
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related spectra. For a — a we considered the QP event 
sample described in the paragraph above. For a-d, we 
have additionally set the deuteron multiplicity to at least 
1. One may see that the a — a correlation function shows 
a narrow peak centered at 78 keV (r=88 keV) which cor- 
responds to the ground state of ®Be ((5=-92 keV) and a 
much broader peak centered at 3.17 MeV (r=4.80 MeV), 
corresponding to the first excited state at 3.03 MeV. For 
a-d we get a well formed peak centered at the excitation 
energy 2.21 MeV (r=0.36 MeV) which definitely corre- 
sponds to the first excited state of ^Li at 2.186 MeV. 
The broadening of correlation peaks is the genuine con- 
sequence of detector finite granularity. The finite granu- 
larity is responsible also for a certain imprecision in the 
determination of the energy. To illustrate this, Fig. [T] 
presents the correlation functions obtained under differ- 
ent hypotheses on the angle under which the detected 
particles were emitted. In the first case (solid symbols) 
we considered that all the particles which hit a certain 
module have been emitted under the angle corresponding 
to the geometrical center of that module. An alternative 
solution is to attribute to each particle a random angle in 
the domains allowed by the geometrical extension of the 
detector (open symbols). As we have noticed that the 
first hypothesis leads to smoother distributions and does 
not introduce extra broadening, the fits are done on the 
corresponding curves. This first hypothesis will be used 
for what follows. Given the complexity of the apparatus, 
one may conclude that the energy calibration is excellent 
and proceed to further spectroscopic analyses. 

To do this, we exploit multi-particle correlation func- 
tions as shown in Fig. [2] for 3-a (a) and 4-a (b) emis- 
sions. In the first case we consider the QP events with 
TOq — 3 and in the second case (poorer statistics) we 
consider QP events with ma > 4. The error bars are cal- 
culated considering only the statistical errors of the cor- 
related spectra. Indeed errors on the uncorrelated spec- 
tra have been reduced to negligible values by increasing 
the number of uncorrelated events as compared to corre- 
lated events. The 3-a correlation function with standard 
event mixing for building the uncorrelated yield spectrum 
(i.e. each particle belongs to a different event, solid sym- 
bols) shows two peaks, at Eex==7M MeV (r=0.18 MeV) 
and 9.62 MeV (r=1.33 MeV). They correspond respec- 
tively to the Hoyle state {E^^p=7.6M MeV, r'=^P=8.5 
eV) and to the complex excited region of ^^C, character- 
ized by the strong EI^p=9M McV, 3" state and by the 
broad £"1^^=10.3 MeV, 0+ state submerging a possible 
2+ state at 9.7 MeV EEll^ ■ To have a first indication on 



the extent to which these states may be considered can- 
didates for a-particle condensation, it is mandatory to 
check whether the decay is simultaneous or if it proceeds 
partially via ^Be. To do this, we adopt the procedure 
proposed in Ref. [2l| and plot with open symbols the 
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correlation function obtained with partial event mixing 
(i.e. the uncorrelated yield spectrum is built by selecting 



FIG. 2: (Color online) Correlation functions in excitation en- 
ergy for 3-a (top) and 4-a (bottom) emission. The distribu- 
tions plotted with solid circles correspond to the case in which 
the uncorrelated yield is built by considering that each par- 
ticle comes from a different event. The distributions plotted 
with open circles correspond to the case in which the uncor- 
related yield is built by considering emission through *Be, 
^^C, 2 X *Be intermediate states. The arrows correspond to 
centroids of Breit-Wigner distributions (solid lines) 



two a particles from the same event and the third one 
from a different event). As one may see, the two peaks 
survive, but while the peak at higher excitation energy 
is almost unchanged, the magnitude of the Hoyle state is 
diminished. This is a qualitative indication that decays 
proceeding via '^Be play a different role in the two cases. 
The bottom panel of Fig. [2] shows that the spectroscopic 
information one may extract out of the 4-a correlation 
functions is much poorer. The first explanation for the 
modest quality of those correlation functions is the re- 
duced statistics generating large error bars. But an even 
more important and fundamental reason is that in this 
energy domain the density of states which decay via a 
emission is high and, not less important, most of them are 
rather broad. Indeed, a numerical simulation where si- 
multaneous 4-a decay is assumed to occur with a branch- 
ing ratio of 1 from each excited state with completely 
missing decay information or which is known to decay 
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FIG. 3: (Color online) Three-a intra-event correlation func- 
tion expressed as a function of average kinetic energy - as^ 
(RMS) of a particles. The uncorrelated yield is built such 
to allow for decay through *Be. The dotted line marks the 
maximum RMS compatible with momentum conservation. 



via a-emission shows that, after filtering by CHIMERA, 
the only state above the 4-a emission threshold which 
may be identified is the one at 15.78 MeV. The solid 
symbols in Fig. [2] correspond to standard event mixing 
for the uncorrelated yield spectrum. Open symbols cor- 
respond to partial event mixing considering decays via 
intermediate ®Be, ^^C and 2 x ^Be. The combined back- 
ground considers the relative probability of these decays. 
In the region of interest, around Eex=i5.1 MeV, a lower 
reduction (as compared to the one of the Hoyle state) is 
observed between the two functions. This seems to be 
an indication that the direct Aa decay is more favourable 
due to larg er decay barriers involved with ^Be and ^^C 
decays [22|. 

To summarize, we demonstrated so far that the 
'^^Ca.+^^C nuclear reaction at 25 MeV/nucleon populates 
excited states of both ^^C and ^^O nuclei which decay by 
3- and, respectively, 4-a emission. Certainly this part 
of information is not sufficient to conclude that those 
states represent a-particle condensation. A major step 
in this direction is to show that the emitted particles 
have practically the same kinetic energy. To do that an 
intra-event correlation method is used. The correlation 
function reads: 



1 + R{aE^,< Ec, >) 



Ycorrjo'Ea, < Eg >) 



(2) 



Here, for a given alpha multiplicity, the numerator is the 
yield of events with given average kinetic energy of alpha 
particles, < Eg >, and given root mean square, (T_e„. 



FiglS] illustrates the intra-event correlation function 
corresponding to 3-a decays from the region of the Hoyle 
state (7.375 MeV^ E^x ^7.975 MeV); the denominator 
is built using partial event mixing. Around the energy 
of the Hoyle state above the 3-a threshold (< E^ >= 
379/3 keV) two peaks are observed. Qualitatively the 
same picture is obtained with standard event mixing and 
also when removing the restriction on the excitation en- 
ergy range. It is remarkable to notice the peak localized 
around < Ea >=110 keV and ctb^ ^ 25 keV. Those 
values correspond, within our energy uncertainties (cal- 
ibration and direction of velocity vectors), to an equal 
sharing of the available energy of the Hoyle state (379 
keV) among the three a particles. The number of events 
which contribute to the peak is 39. So we can consider 
that at least this number of condensate states have been 
populated in the present experiment. The peak localized 
around < Ea >=130 keV and around 90 keV cor- 
responds to the sharing of the available energy between 
the two as of ^Be and the remaining a of 191 keV. The 
number of such events is 85 and must be also considered 
as a lower limit. Indeed, depending on their velocity and 
polar angle in the laboratory, ^Be from the Hoyle state 
can decay by firing a single telescope. For the sake of 
completeness, we mention that the total number of events 
with (m8_Be=l, iTia = 1) in the region of the Hoyle state is 
900 while that corresponding to TOq, = 3 is 1072. Numer- 
ical simulations, filtered by the multi-detector replica, 
confirm that the observed peaks in fig. [3] correspond to 
3a decay (RMS<20 keV) and ®Be-a decay, with the two 
a flying along the initial *Be direction (< Ea >~ 130 
keV, RMS ~85 keV). The other extreme case, *Be-a de- 
cay, with the two a flying perpendicularly to the initial 
^Be direction would fall at (< Ea >^ 90 keV, RMS -60 
keV). The branching ratio between the two decay pro- 
cesses, very delicate to estimate, can not yet be given. 
Concerning the 3-a decays of the broad 9.64 MeV state, 
we mention that no peak at low ctb^ was observed in the 
intra-event correlation function. The much poorer statis- 
tics we have for 4-a does not allow a similar analysis. 
Nevertheless, we mention that there are four events com- 
patible with the above stated criteria and corresponding 
to the excited state at 15.1 MeV of ^^O [isf . 

In conclusion, the nuclear reaction ^°Ca-|-^^C at 25 
MeV/nucleon bombarding energy was used to produce 
states theoretically predicted as a-particle condensate 
states. Supposing that equal values of kinetic energy of 
the emitted a-particles represent a sufficient criterion for 
deciding in favor of a-particle condensation, we found 39 
events corresponding to the direct de-excitation of the 
Hoyle state of ^^C which satisfy the criterion. It also 
clearly appears that the major part of Hoyle states deex- 
cite via the ^Be-l-a channel. Though the poor statistics 
prevents drawing a definite conclusion in the case of ^^O, 
we mention the existence of a small number of events 
originating from the decay of the 15.1 MeV state which 



fulfill the present criteria. An experiment with higher 
statistics is planned to better study the ^^0 case. 
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